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Abstract 

The difference of temperatures between an impurity and the surrounding gas in an open vibrated 
granular system is studied. It is shown that, in spite of the high inhomogeneity of the state, the 
temperature ratio remains constant in the bulk of the system. The lack of energy equipartition is 
associated to the change of sign of the pressure diffusion coefficient for the impurity at certain values 
of the parameters of the system, leading to a segregation criterium. The theoretical predictions 
are consistent with previous experimental results, and also in agreement with molecular dynamics 
simulation results reported in this paper. 
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Granular systems, i.e. assemblies of macroscopic grains dissipating part of their kinetic 
energy during collisions, exhibit a rich phenomenology with many qualitative differences as 

n 

compared with molecular systems 111. Among these are the breakdown of energy equipar- 

u n . n n 

tition |2J, |3J and particle segregation in agitated mixtures j^, |5|| . The aim of this work is to 
investigate both effects for a dilute mixture in the tracer limit, i.e. very low mole fraction 
of one of the species, showing that there is a close relationship between them. 

Different species of a granular mixture have different granular temperatures, defined from 
the average kinetic energy of the particles. A prediction for their ratio in a dilute binary 
mixture of smooth inelastic hard spheres has been derived from the inelastic Boltzmann 
equations describing the time evolution of the distribution functions of the species (J, and 
its accuracy has been verified by Molecular Dynamics simulations for the particular case 
of an isolated homogeneous system 7]. In this Letter, it will be shown that kinetic theory 
also accurately predicts the temperature ratio in the case of a highly inhomogeneous driven 
system. This ratio turns out to be constant in the bulk of the system, consistently with 
the experimental results reported in ^j]. As in the homogeneous state, deviations from 
equipartition depend on the mechanical differences between the species and the degrees of 
inelasticity of collisions. 

Particle segregation is the demix of a granular mixture when shaken. Usually, the larger 
particles are observed to rise (Brazil-nut effect), although under certain conditions they 
can also tend to descend (reverse Brazil-nut effect). Several mechanisms, corresponding to 
different scenarios, have been proposed to explain both behaviors jj, 0, 0], although the 
phenomenon is far from being fully understood. Here, segregation will be investigated in 
he context of hydrodynamics for a dilute granular mixture as derived from kinetic theory 
lp| . The relative position of the tracer component with respect to the excess component 
is determined by the sign of the pressure diffusion coefficient. While in a molecular gas 
this sign is fixed by the mass ratio of the particles of the components, for a granular gas it 
also depends on the temperature ratio. Due to the lack of energy equipartition, if follows 
that the criterion for segregation is rather complicated, involving all the parameters of the 
mixture. 

The model system considered is a low density gas of smooth inelastic hard spheres (d = 3) 
or disks (d = 2) of mass m and diameter a, and one impurity of mass mo and diameter cr - 
This is formally equivalent to the tracer limit for the impurity component. Inelasticity of 
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collisions is specified by two independent constant coefficients of normal restitution, a and 
QiO) referring to gas-gas and impurity-gas collisions, respectively. The system is in presence 
of a uniform external field of the gravitational type, so each particle is submitted to a force 
per unit of mass given by —goS z , where go is a positive constant and e z a unit vector in the 
direction of the z axis. Energy is continuously supplied to the system through the bottom 
wall located at z = that is vibrating with small amplitude and high frequency. There is 
no upper wall, i.e. the system is open. 

Under the above conditions, the system exhibits an inhomogeneous steady state with 
gradients only in the z direction and vanishing velocity field. In the case of a one-component 
system, it has been verified that the hyd ro dy namic profiles away from the walls are well 



described by the hydrodynamic equations 



roclynai 

bills. 



12| . It is assumed that these profiles are not 



affected by the inclusion of the impurity. Consider the local temperatures of the gas and 
the impurity, T(z) and T (z), defined in the usual way from the respective mean square 
velocities (with the Boltzmann constant set equal to unity). A formal relation between 
both temperatures can be derived if the existence of a hydrodynamic regime is assumed 
and, moreover, that the corresponding "normal" solution of the Boltzmann equation for 
the impurity can be generated by an extension of the Chapman-Enskog method jloL [l^ . 
Then, it is obtained that the lowest order in the gradients of the cooling rates for the gas, 
£(°\z), and the impurity, (o°\z), must be equal. These rates are nonlinear functionals of 
the zeroth order distribution functions, and can be estimated at good approximation using 
Maxwellians, with the results [3| 
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where n(z) is the number density of the gas, v g (z) = [2T{z)/m] 1 / 2 , h = m(l + a )/2(m+m ), 
Uq is a dimensionless impurity-particle collision rate, 

8h7T (d-i)/2 /w ^d-i 



u T(d/2)d \<t. 
with a = (cr + o"o)/2, and <fi the ratio of mean square velocities, 

mT (z) 
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Equating Eqs. (JTJ and (J2J) provides the equation for 



This gives a cubic equation which has a unique real, positive solution for all the allowed values 
of h and /3. For elastic collisions one gets <p — w/mp, as required by energy equipartition. 
Moreover, in the limit of an elastic gas, a — 1, but inelastic collisions between the intruder 
and the fluid, aio < 1, the expression Obtained i n re , H is _d. Toe oe hOT of 
the solution in the limit <fi — > in the context of the isolated homogeneous state has been 
analyzed in fl6j ]. where it has been shown that a qualitative change similar to a second order 
phase transition occurs, but this will not be discussed here. 

From Eq. (JHJ) it follows that (f> does not depend on z. Therefore, the temperature ratio 
T (z)/T(z) not only differs in general from unity, but remains constant along the system. 
Moreover, it does not depend either on the precise way in which the wall is being vibrated, 
as long as the assumed steady state is reached. Of course, the result only holds in the 
bulk region of the system, where a hydrodynamic description is expected to apply. These 
predictions are in qualitative agreement with the experimental results found by Feitosa 
and Menon |3j, although it must be realized that their experiments are not carried out in 
the tracer limit. On the other hand, these authors did not find any dependence of the 
temperature ratio on the inelasticity of the grains, while Eq. (jSJ) predicts a rather strong 
one. 

To check the above results, we have performed molecular dynamics (MD) simulations 
of a system of N inelastic hard disks (d = 2). The wall at the bottom vibrates with a 
sawtooth profile with velocity v w and collisions of the particles with it are elastic. Moreover, 
the amplitude of the wall motion is much smaller than the mean free path of the particles 
next to it, so in practice the position of the wall can be taken as fixed at z — 0. Periodic 
boundary conditions are used in the direction perpendicular to the field. The units are 
defined by m = o = g = 1 . In all the results to be reported in the following, the parameters 
determining the state of the gas have been kept constant, namely N = 359, a = 0.95, 
v w = 5, and a width S = 50. For these values, it is verified that the system is fluidized and 
actually reaches the steady state assumed above. In Fig. [H the temperatures profiles of the 
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FIG. 1: Dimensionless temperature profiles of the gas (solid line) and the impurity (symbols) for 
mo/m = 1/2. The different symbols correspond to different values of the restitution coefficient 
Q!q ; as indicated. The values of the parameter h, defined in the main text, are also indicated for 
reference. 

impurity for mo/m = 1/2, do = 1, and several values of «o are shown. Also, the temperature 
profile of the gas (solid line) is included. It is clearly observed that energy equipartition is 
not verified, and that the temperature difference depends on a - The ratio between the 
impurity and gas temperatures for the same states is plotted in Fig. |21 In agreement with 
the theoretical prediction, the ratio remains constant, aside from statistical fluctuations, in 
the interior of the system, even in the region where each of the partial temperatures presents 
large gradients. In fact, the kinetic boundary layer next to the vibrating wall is quite narrow. 
Similar behaviors have been found for other values of m , namely m /m = 0.75, 1 and 2. 

The quantitative comparison of the simulation results with Eq. (j5J) is presented in Fig. 
EJ where <fi is plotted against the parameter h. The solid line is the solution of the equation 
and the different symbols correspond to different values of the mass ratio as indicated. A 
very good agreement is obtained. It is worth to stress that the theory remains accurate up 
to very small values of the coefficient of restitution ato, indicating the validity of the kinetic 
and hydrodynamic descriptions even for very strong dissipation. Also, for later discussion 
note that <fi can be larger than unity even if m/m < 1, and viceversa. 
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FIG. 2: Profiles of the ratio between the impurity and the gas temperature for the same values of 
the parameters as in Fig. ^ 



Consider next the number density for the impurity in the steady state, n (z). It is 
determined by the condition that the associated flux j z must vanish. To first order in the 
gradients, it is 

77777 77? 

jz = -m Dd z x - 7^D>d z T - —D p d zP (z), (7) 

where Xq = n^/n, p = nT is the pressure, D the diffusion coefficient, D' the thermal 
diffusion coefficient, and D p the pressure diffusion coefficient. Explicit expressions for these 
transport coefficient can be derived by the Chapman-Enskog procedure in the first Sonine 
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where A = mj [m + mo) and v e is the elastic limit, 
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FIG. 3: Ratio of the mean square velocities as a function of the parameter h defined in the main 
text. The solid line is the theoretical prediction given by Eq. © and the symbols are from the 
simulations for different values of the mass ratio, as indicated . 

Although the equation j z = can be (numerically) integrated to get n Q (z), given that 
the hydrodynamic profiles for the gas are known, we will restrict ourselves to a qualitative 
property of n (z). We will determine the position of the maximum of Uq{z) relative to that 
of n(z). That means that we consider situations in which the latter exists, which is always 
the case for large enough number of particles N Let us denote by z m the height at 

which this maximum occurs. It has been established that at t his p oint the temperature of 
the gas is a decreasing function of z, i.e. dT/dz < at z = z m [llj. Then, it is obtained 
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If the right hand side of the above equation is positive, uq(z) is still growing at z — z m , 
following that its maximum occurs at 2o,m > z m . In the opposite case, the maximum of Uq{z) 
takes place at z ,m < z m . Note that this sign is the same as that of the pressure diffusion 
coefficient D p . The general discussion is rather complicated, given the large number of 
parameters involved [13|. For the sake of simplicity, we are going to consider the region of 
the parameter space verifying (3 < h/2, which includes all the MD simulations reported here. 
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FIG. 4: Normalized height distributions for the gas (solid line) and the impurity (symbols) for the 
same systems as in Fig. ^ 

For this range of values, it is easily verified that 1 — /2u + /2u 2 > 0, and the sign 
of dlnno/dz at z = z m is determined by the value of 0. For <fi > 1 (0 < 1) the derivative 
is positive (negative) and the position of the density maximum for the impurity is higher 
(lower) than that for the gas. As mentioned above, the value of <ft can be quite different 
from the mass ratio due to the different temperatures of the impurity and the gas. Since 
the value of is a function of all the parameters of the system, the segregation criterion, 
understood as a rule determining the relative position of the impurity with respect to the 
gas, also involves all those parameters. 

The MD density profiles for the same systems as in Fig. ^ are given in Fig. 0] Again, the 
solid line is the density profile of the gas, normalized to unity for comparison purposes. As 
the restitution coefficient a.Q increases, the maximum of the density profile for the impurity 
moves towards higher regions and the profile becomes wider, showing the tendency of the 
impurity to rise. The same effect is observed for the other values of the mass ratio considered. 
Since the value of increases as «o increases keeping the remaining parameters constant, 
this behavior is in qualitative agreement with the above theoretical predictions. 

The distributions for the position of the impurity become quite flat as the value of a 
increases, rendering quite difficult to identify their maximum from the MD results. Then, 
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FIG. 5: Ratio between the impurity and the gas center of mass positions versus the ratio of the 
mean square velocities 4>. The different symbols correspond to different values of the mass ratio as 
indicated. 

what has been measured is the center of mass position for both the gas, z cm , and the 
impurity, z . In Fig. El the ratio z Q /z cm versus <fi is shown. It is seen that the segregation 
criterium derived above is fairly obeyed, in spite of the fact that the shape of the density 
distribution clearly indicates that the position of the maxima are rather different from those 
of the centers of mass. Similar results have been found for other values of the parameters, 
including different sizes of the impurity jl^ ]. 

To conclude, it must be stressed that the present work deals with a dilute granular 
mixture in the tracer limit. At higher densities, other segregation mechanisms, such as 
those discussed in j^] and become important. On the hand, there is no reason to expect 
that the one identified here fails to be relevant at those densities. 
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